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We present a construction scheme for quantum error cor- 
recting codes. The basic ingredients are a graph and a fi- 
nite abelian group, from which the code can explicitly be ob- 
tained. We prove necessary and sufficient conditions for the 
graph such that the resulting code corrects a certain number 
of errors. This allows a simple verification of the 1-error cor- 
recting property of fivefold codes in any dimension. As new 
examples we construct a large class of codes saturating the 
singleton bound, as well as a tenfold code detecting 3 errors. 



I. INTRODUCTION 

From the beginning of Quantum Information Theory 
it was recognized that error correcting codes play a cru- 
cial role. On the one hand it was clear that without 
error correction decoherence effects could easily annihi- 
late the gain in computing time promised by the new fast 
quantum algorithms. On the other hand, the no-cloning 
theorem seemed to forbid at least the most naive ap- 
proach to classical error correction for noisy channels, 
e.g., sending each bit three times and taking a majority 
vote at the output of the channel. Clearly, this simple 
scheme reduces classical errors with small probability of 
order e to order e 2 . The cloning required for sending 
" the same bit" three times rules out direct quantum ana- 
logues of this scheme. It was therefore an important step 
to realize that quantum mechanics had other, more 
subtle ways of "distributing" quantum information over 
several channels to stabilize against errors. One problem 
with the known schemes of quantum error correction (e.g. 
HQ), however, is that they tend to be subtle indeed, and 
the verification of their error correcting capabilities often 
requires a lengthy computation. It is therefore desirable 
to find new, perhaps simpler ways of constructing error 
correcting codes, on which more direct intuitions might 
be built. 

In this paper we propose a new scheme for construct- 
ing quantum error correcting codes, which has some of 
these advantages. The ingredients of our construction 
are a graph and a finite abelian group. The order of 
the group determines the type of systems for which er- 
rors are corrected so that, e.g., the two-element group 
corresponds to the qubit case (Compare [^|-|| for other 
constructions of non-binary codes). The number of ver- 
tices of the graph determines how many input systems are 
coded into how many output systems. From the edges of 
the graph one can then read off an explicit expression 
for the code. However, not every graph corresponds to a 



good code, and we will discuss the condition for the code 
to correct a certain number of errors. In the simplest 
case, the fivefold code P-|TT[| (for qubits as well as higher 
dimensional systems), it can be verified in a few lines that 
any two errors are corrected. We also give an example of 
a more complex tenfold code detecting 3 errors. 

As we are going to discuss in a following paper in 
more detail, it turns out that the codes which can be 
achieved by our method are stabilizer codes. There are 
various efficient methods for constructing stabilizer codes 
P,0JT3-[l3]. However, we think that, compared to previ- 
ous stabilizer constructions, our technique has some in- 
teresting new features. 

• Often the condition for error correction can be 
proved for many groups simultaneously, so that one 
gets code families for systems of variable sizes. 

• The geometric intuitions about graphs may become 
helpful for finding new constructions. 

• Our codes have the property that in their natu- 
ral basis all matrix elements of the coding operator 
have the same modulus (Hadamard form). This is 
prima facie opposite to the usual goal of getting as 
few non-zero matrix elements as possible. However, 
the latter can be achieved for our codes by a dis- 
crete Fourier transform applied to some channels. 
Moreover, the Hadamard form appears to be an in- 
teresting normal form for the codes which can be 
written in both ways. 

• For some codes it is possible to exchange some in- 
put vertices with some output vertices, while re- 
taining the error correction property. This kind of 
symmetry is much harder to see in the usual sta- 
bilizer constructions, and may prove to be helpful 
in coding problems with additional inputs and out- 
puts, such as the internal state of the codig device 
in convolutional codes. 

The paper is organized as follows: We begin by de- 
scribing the general construction of the coding operator 
within Section D. In Section III we recapitulate the Knill- 
Laflamme condition for error correction and adapt it to 
our particular type of codes, resulting in a necessary and 
sufficient condition for a graph to generate a quantum 
error detecting code. The remaining sections contain ex- 
amples of codes constructed in this way. In Section |^ 
we show that it becomes simple indeed to verify the five- 
fold quantum codes. In Section ^ we demonstrate that 
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for a given number e of errors and number k of inputs, 
there is a graph generating an infinite code family using 
4e + k output systems, i.e., a family of codes saturating 
the singleton bound. Finally, in Section VI we construct 
a code with 1 input and 10 outputs, detecting 3 errors 
for arbitrary system size. 



II. BASIC CONSTRUCTION 

Every code we construct is completely determined by 
the follow ingredients: 

• An undirected graph T with two kinds of vertices. 
We distinguish the set X of input vertices and the 
set Y of output vertices. The links of the graph are 
given by the coincidence matrix of the graph, which 
we will denote by T for short. Its matrix element 
r(zi,z 2 ) is 1 iff the vertices zi,z 2 £ (X U Y) are 
linked, and otherwise. More generally, we allow 
weighted graphs, whose incidence matrices have ar- 
bitrary integer entries, apart from the constraints 
r(zi, z 2 ) = T(z 2 , zi) and F(z, z) = 0. 

• A finite abelian group G with a non-degenerate 
symmetric bicharacter. 

By definition, a bicharacter is a function x:GxG->C 
such that x(g + h,g') = x{9i9')x{h,g') and a similar 
condition holds for the second argument, which is also 
implied by the assumed symmetry x(ffjS') = x(tfi9)- 
We also assume non-degeneracy in the sense that 
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\G\ S(g') 



\G\ for 5 ' = 
for g> ^ 



(1) 



Note that since every g £ G has finite order, \{9i s') is 
always a root of unity, and x(<7, <?') = x{~9i </)■ For G — 
Z p , the cyclic group of order p, the standard bicharacter 
is given by 



X(g,h) = exp(— ghj 



(2) 



where g, h are integers representing their class modulo 
p. Since every finite abelian group is a direct product of 
cyclic groups, this also shows the existence of bicharacters 
for any such group. 

The input and output systems of the code are labeled 
by X and Y. They are all of the same type, i.e., they 
are described by the same Hilbert space 7i — L 2 (G). 
This is the space of all functions tp : G — > C with scalar 
product {(pfijj) = | C | 1 J2 g 4>{g)^p{g). For compactness 
of notation we write such normalized sums as integrals. 
Hence the scalar product becomes J dg cf>(g)ip(g). The 
combined input system is thus described in the |A|-fold 
tensor product TL® X = L 2 (G X ). Vectors in this space 
are functions of \X\ variables, one variable g z for every 



z £ X. The entire collection of variables will be denoted 
by g x . The error correcting code will be an isometry 



v r : L 2 {G X ) -» L 2 {G Y ) 



(3) 
(4) 



where vr under the integral denotes the integral kernel of 
the operator vr- This kernel depends on both input and 
output variables, which are combined into a collection of 
variables g z , one for each vertex z £ X U Y of the graph. 
The core of our construction is an explicit expression for 
this integral kernel: apart from an overall normalization 
factor, it will simply be a product of phases, with each 
factor corresponding to a link of the graph: 



v r [«^] = \G\^ 2 J] X{9^f {z ' z ' 

{z,z>} 



(•5) 



where the product is over all 2 element subsets {z, z'} C 
(X U Y). Thus for an ordinary graph {T(z,z') = 0,1), 
this is the product of all x(ffz,<?z')j f° r which z and z 1 
are linked. The remarkable property of such codes is 
that apart from the normalization factor the kernel is 
everywhere of modulus 1. When G = Z p is cyclic, and x 
is given by Equation (||) , we can write the phase in more 
compact form as 



X\JY} 



= \G\\ X \' 2 



exp 



m „X\JY 



f 7TC 

V p 



(J 



XUY 



(0) 



where "•" denotes the product of integer valued matrices 
and vectors. Note that every term in the sum g XuY ■ T ■ 
qXuY occurs twice, which we compensated by a factor 
1/2- ^ 

This completes the construction of the operator v from 
the defining ingredients listed at the beginning of this 
section. Of course, in general this will not be an error 
correcting code nor even an isometry. The conditions for 
this will be studied in the following section. 



III. THE CONDITION FOR ERROR 
CORRECTION 

A general characterization of quantum error-correcting 
codes has first been worked out by E. Knill and R. 
Laflamme [|L7| . We briefly review here the main aspects, 
and adapt the condition the particular case of codes con- 
structed as in the previous section. In this theory a quan- 
tum code is an isometry v : TL — > IC from the "input 
Hilbert space" TL to the "output Hilbert space" JC. Thus 
an input density operator is transformed by coding into 
vpv*, which is a density operator on IC. The output of 
the coding is then passed through a noisy channel. The 
noise is described by a certain class of errors, which are 
represented by a linear subspace £ of operators on IC. 
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The channel is thus represented by a completely positive 
linear map of the form 



T{p) = Y,F aP F* 



(7) 



where F a G £. and are chosen such that the output is 
always normalized. The isometry v is said to be an error 
correcting code for E , if there is a completely positive 
"recovery operator" R such that 



R(T(vpv*)) =p 



(8) 



for all density operators on TC. By the theory of Knill 
and Laflamme JT^I this is equivalent to the factorization 
condition 



np x ,F*F^ 2 )=u>{F*Fp) (^2} 



(9) 



where uj(F a , Fp) is a factor independent of the arbitrary 
vectors tj)x,ip2- As in much of the literature on codes we 
will consider here a specific type of errors, namely errors 
happening only on a small number of outputs of the code. 
Thus the tensor product structure H® Y — L 2 {G Y ) of the 
output space becomes important. Let 2l(-E) denote the 
set of all operators on L 2 (G Y ), which are localized in 
E C Y, i.e., which are the tensor product of an arbitrary 
operator on H® E with the identity on 7i® Y \ E . We say 
that a code corrects e errors, if F a , Fp in (^) may be 
chosen arbitrarily in the linear span of U|B|< e Note 
that the operators F*Fp appearing in the scalar product 
@ can then be localized on arbitrary sets of 2e elements 
and any operator with such localization may be written 
as a linear combination of such F*Fp. It is therefore 
convenient to introduce the following terminology: we 
say that the code v detects the error configuration E C V, 
if 



(vipi,Fvip 2 ) =oj(F) (ipuipz) , 



(10) 



for all F £ 2l(-E). Then a code corrects e errors, iff it 
detects all error configurations E C Y with \E\ < 2e. 

We will now adapt these conditions to operators vr of 
the special form (||) . Consider a fixed error configuration 
E CY, and let I = Y \ E. Then if F is an operator F on 
TC® E , with integral kernel F[g E , h E \, the integral kernel 
of v*Fv is 

vlFvr\g x ,h x ] = J dg E dg I dh E v r [ 5 * 5 B , 5 '] 

X^fcV^,)!^']. (11) 

This must be a multiple of the identity for every choice 
of F. Choosing, in particular, a rank one operator 
F = \g E )(h E I we find that error detection for the configu- 
ration E is equivalent to the property that the correlation 
function 



^r, E] [9 XUE ^ XuE } := J dg 1 v r [ 5 * 

factorizes in the following manner 

^rM9 XUE ,h XuE ] = C(g E ,h E ) 5{g x h x ) , (13) 

where S(g x ) is defined to be 1 if g x = for all x G X, and 
zero otherwise, and C(g E , h E ) is a factor independent of 
the input variables g x , h x . 

For two subsets K, L of X U Y we denote by the 
group homomorphism from G L to G K which can be de- 
rived from the corresponding submatrix of the incidence 
matrix T by the prescription 



T l9 L ■= (E r ( fc ' Z )3/ 



leL 



(14) 



and the following condition is necessary and sufficient for 
quantum error detection: 

Theorem III.l Given a finite abelian group G and a 
weighted graph T as in the basic construction. Then an 
error configuration E C Y is detected by the quantum 
code vr iff the system of equations 



1 xuE a — u 



with I = Y \ E implies that 



d x = and T x d E = 



(15) 



(16) 



The proof of Theorem III.l is given in Appendix ^ 
and, in view of this result, we discuss some examples. 



IV. EXAMPLE: THE FIVEFOLD WAY 

The first example of an optimal quantum error cor- 
recting code correcting all one bit errors was the famous 
five qubit code by |l(|. The original code is not easy to 
verify, so it is gratifying to see that our construction pro- 
duces such a code, which can be verified in a few lines. 
Moreover, our construction works simultaneously for all 
groups G, and is hence not restricted to qubits. Five- 
fold codes for higher dimensional systems have been con- 
structed before jn]], and if we believe a recent result by 
Rains fl8|| , the qubit code is essentially unique anyway. 
Hence this section has mainly illustrative character. 

Consider the graph in Figure ^ where the central vertex 
"0" is the input vertex, and the remaining five are the 
output vertices. 
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FIG. 1. The left figure shows a graph for a fivefold code, 
where the central vertex is the input vertex. The two figures 
on the right side are the relevant 2-error configurations. 



We will verify the condition of Theorem [II. 1 in a par- 
ticularly strong form. Namely we will show that, for 
every 2-element error configuration E, 



r I JXUE 
1 XUE a 







implies 



jXuE 







(17) 



The error configuration is a 2-element subset of the out- 
put vertices 1, ...,4, and for the purpose of verifying crite- 
rion ([l7]) the input vertex plays exactly the same role as 
an error. It is clear by symmetry that only the two con- 
figurations for XUE shown by filled dots in Figure [j] need 
to be considered. Now the condition Y I XUE d XUE = is a 
set of equations, one for each "integration vertex" y E I: 
For each vertex y we have to sum the d x for all vertices 
of x G X U E connected to y, and equate to zero. (In a 
weighted graph, we would have to sum with coefficients 
given by the matrix V). The following is a table of equa- 
tions arising in this way for the first error configuration, 
XI) E = {0,1,2}: 



Vertex y 


Equation 


3 
4 
5 


d + d 2 = 
do = 
d + di = 



Clearly, this implies do = d\ = d% = in any abelian 
group. Similarly, for the second error configuration X U 
E = {0, 1, 3}we get the equations 



Vertex y 


Equation 


2 


do - 


h di + d 3 = 


4 




d + d 3 = 


5 




d + di = 



which once again implies do = d\ = d 2 = 0. This con- 
cludes the verification that the code associated with the 
graph in Figure |l|, and an arbitrary finite abelian group 
G, detects any two errors, hence corrects one error. 

In fact, we proved a little bit more than that. The 
essential part of the proof was to look at certain 3x3- 
submatrices of the 6 x 6-matrix T, namely those corre- 
sponding to an off-diagonal block in the partition of the 
vertices {0,1,..., 5} into two disjoint subsets X li E and 
/, and to show that each such submatrix is non-singular. 
Regarded in this way, it becomes irrelevant to which of 
the two sets in the partition the input vertex "0" be- 
longs, so we showed that any vertex, even a peripheral 
one, may be taken as input vertex, and we still get a 
1-error correcting code. 

This may seem like a rather strong property of the 
graph we chose. However, there is (exactly) one other 
graph with six vertices, which produces in the same way 
a 1-error correcting code for arbitrary choice of input 
vertex and abelian group G. This is shown in Figure^. 





FIG. 2. An alternative graph for a fivefold code which cor- 
rects one error. 



V. CODES SATURATING THE SINGLETON 
BOUND 

In this Section we briefly discuss one natural general- 
ization of the idea emerging in the previous section. For 
definiteness, let us fix G as a cyclic group of prime order 
d, so Z,d is a field, and choose an integer m. We then 
ask for symmetric 2m x 2m-matrices Y with integer en- 
tries (or, equivalently, entries in Z<j) with the following 
property: for any m-element subset / C {1, . . . , 2m} the 
m x m-submatrix Ty with i € J, j ' ^ I is invertible in the 
field Id- For the purpose of this section, let us call such 
a matrix strongly error correcting for the prime number 
d. 

What codes can we get from such matrices? Just as in 
the previous section, let us specify any set X of k < m 
vertices as input vertices, and call the remaining n — 
(2m— k) ones output vertices. Then for any configuration 
of e = m — k errors, the set XUE and its complement I 
will have exactly m elements. By assumption, the strong 
form of the error correcting condition fll7j ) is satisfied, 
hence the code detects e errors. These parameters satisfy 
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n = 2e + k 



(18) 



i.e., the general inequality n < 2e + k, known as the 
singleton bound Jl7| is satisfied with equality for any such 
code. 

How can one get strongly error correcting matrices T 
in a practical way? Here is a procedure, we found easy to 
work with for small to, using a symbolic algebra program. 
First, we introduce variables for each matrix element 
with i > j, and compute the determinants of all off- 
diagonal to x TO-submatrices as symbolic expressions in 
these variables. As we go along fixing integer values for 
these Tij, the determinant expressions become simpler, 
and in some cases factorize. Each of these factors has 
to be kept non-zero by the next choice of a -value. 
Finally, we end up with an integer matrix, whose off- 
diagonal to x TO-submatrices all have non-zero integer 
determinant. Then, for any prime d, which does not 
divide any of these integers, we have solved the problem. 

It is natural to begin this process by setting as many 
weights as possible equal to zero. It is easy to see that 
Theorem [II. 1 does not allow too many r(fc, I) = be- 



is just the number of lines meeting a particular vertex as 
one can see from Figure |[ 



cause an entire row of zeros in the matrix T I XUE leaves 
one of the difference variables unconstrained. Similarly, 
in the condition for strong error correction it is clear that 
no off-diagonal submatrix should have a row of zeros, i.e., 
each one of the 2m vertices must be connected to at least 
to other vertices. The graph in Figure || is as sparse as 
possible under these constraints (to = 4), and was the 
starting point for a search for non-zero weights, as de- 
scribed above, resulting in the matrix 



r = 



/0 


1 





1 


1 


1 


o \ 








1 


1 


1 





1 


1 








2 





-1 


1 


1 











1 


2 


-2 


1 1 


2 











-2 





1 1 





1 











-1 


1 


-1 


2 


-2 











V° 1 


1 


-2 





-1 





o / 



(19) 



This matrix can either be used to get codes detecting 
three errors on an arbitrarily chosen single input vertex, 
or as a code detecting two errors (or correcting 1) on 
two arbitrarily chosen inputs. The set of determinants is 
{-11, -8, -5, -4, -2, -1, 1, 2, 4, 5, 8, 9}, so this will work 
for any prime d not in in the set {2, 3, 5, 11}. By fixing the 
choice of the input vertices, we may restrict to a smaller 
set of partitions (the input vertices always belong to the 
same set), hence we get fewer constraints. For example, 
the code with input vertices {1,2} has no relevant sub- 
determinant containing a factor 3, so the resulting code 
corrects one error on arbitrary pairs of 3- level systems. 

Within the above example, the number of matrix ele- 
ments set to zero is maximal, namely v% = 4 for each row 
i = 1, ■ • • , 8. Looking at the corresponding graph, Vi = 4 




FIG. 3. Graphical representation of the weighted graph in 
Equation (19). Edges without label have weight 1. 

Strongly error correcting matrices exist in any dimen- 
sion, so the code parameters saturating the singleton 
bound ( |l8| ) can be chosen arbitrarily, if the dimension 
d of the one-site system is taken to avoid a certain fi- 
nite set of primes (Compare also ||). The argument is 
quite simple: consider the to x to- subdeterminants of 
symmetric 2m x 2TO-matrices as a family of polynomi- 
als f a , a — 1, • • ■ , i 2 ™) in i 2 ™) variables. None of these 
vanishes identically and since Z is an integral domain, 
which in contrast is wrong for finite fields, the product 
polynomial Y[ a fa ^ is non zero jl^, p. 106]. This 
implies that there exists an integer tuple T of arguments 
that ria/atT) 7^ 0. Thus we have shown the following 
statement: 

Proposition V.l For each number e of errors, there 
exists a prime d and a weighted graph T such that the 
quantum code, associated with the weighted graph T, is 
a quantum error correcting code, which encodes k d-level 
systems into 4e + k p-level systems, and which corrects e 
errors. 



VI. A 10-FOLD QUANTUM ERROR-DETECTING 
CODE 

In this section we present more complex example for 
a graph, which yields for every finite abelian group G a 
10-bit code, detecting 3 errors is given by Figure ||. At 
the first look, this graph looks rather complicated, but it 
can be described in a simple fashion, by looking at the 
graph for the fivefold code in Figure [j]. 




FIG. 4. The graph for a tenfold code. 
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Namely, the graph, given by Figure [|, can be obtained 
as follows: The output vertices {1, 2, 3, 4, 5} of the graph 
in Figure [l] are replaced by pairs: 1 1— > {1, 2}, 2 i— > {3, 4}, 
. . . , 5 i— > {9, 10}. Each output vertex is connected with 
the following vertices: The central input vertex 0, the 
vertex belonging the same pair, and all output vertices 
belonging to neighbored pairs. 

The symmetry of this graph can efficiently be used to 
check that each error configuration with 3 errors can be 
detected. As it is depicted by Figure |L there are only 
four 3-error configurations to distinguish: 

• All errors occur within different pairs and all these 
pairs are neighbored (first graph in Figure ^|) . 

• All errors occur within different pairs and only two 
of these pairs are neighbored (second graph in Fig- 
ure ||). 

• One pair is totally affected by errors and the re- 
maining error occurs within a pair which is not 
neighbored (third graph in Figure ^) . 

• One pair is totally affected by errors and the re- 
maining error occurs within a neighbored pair (4th 
graph in Figure |5|) . 




Proposition VI. 1 For each finite abelian group G, the 
quantum code, which is associated with the graph given 
by Figure is a quantum error- detecting code, encoding 
1 input system into 10 output systems and detecting 3 
errors. 

Proof: Suppose that each error occurs in different pairs 
and all these pairs are neighbored (first graph in Figure 
|5|), e.g. the error configuration {1,3,5}. Then we pro- 
ceed as in the example for the fivefold code to obtain the 
system of equation (|l5|): 



Vertices y 


Equation 


2 
4 
6 

7 and 8 
9 and 10 


d + di + d 3 = 
d + di + d 3 + d 5 = 

do + G? 3 + d 5 = 
do + d 5 = 
d + di = 



and we conclude that do — and ^1+^3 + ^5 = are 
true. Thus the equations ( fHf ) are fulfilled and the corre- 
sponding error configuration is detected. 

Analogously, one checks that each error configuration, 
where each error occurs in a different pair and only two 
of these pairs are neighbored (second graph in Figure 
is also detected. 

We now consider an error configuration, where one pair 
{1, 2} is totally corrupted and the remaining error vertex, 
e.g. 5, is contained within a non-neighbored pair (third 
graph in Figure ^|). We obtain for the system of equations 





Vertices y 


Equation 


3 and 4 
6, 7 and 8 
9 and 10 

7 and 8 


d Q + di + d 2 + d 5 = 
do + d 5 = 
d o + d 1 +d 2 =0 
do + G?5 = 



and we conclude again that do — and d\ + d 2 + ds = is 



true. The equations (16) arc fulfilled which implies that 
the corresponding error configuration is also detected. 

Finally, it is a bit more straightforward as in the pre- 
vious case to show that an error configuration, where 
two errors occur within one pair and the remaining er- 
ror within a neighbored pair (4th graph in Figure g|) is 
detected. □ 



VII. CONCLUSION AND OUTLOOK 



FIG. 5. A complete family of inequivalent error configura- 
tions for the graph in Fig. 3. 



In this paper we had to limit the exploration of our 
construction scheme to a few examples. A more system- 
atic investigation is, of course, under way. Some of the 
issues in this investigation are the following. 
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We already mentioned that all codes constructed 
from graphs are stabilizer codes. This is verified by 
explicitly constructing a group of unitaries, com- 
posed out of shifts and multiplication by characters, 
leaving all vectors in the code space invariant. The 
converse of this statement is not so clear, i.e., how 
to embed the usual stabilizer code constructions 
into our scheme, and to characterize the subset of 
codes for which this is possible. 



• We have seen in Sec. IV that different graphs gener- 
ate 5-qubit codes, although such codes are presum- 
ably unique up to local transformations. It would 
be helpful to characterize the local unitary trans- 
formations taking one graph code into another, and 
to study the relationships between the resulting 
graphs. The Rains invariants f20|-f22| for graph 
codes can be computed relatively easily, and should 
help to decide such isomorphism issues. 

• From Sec. fv| it is clear that the singleton bound be- 
comes easier and easier to satisfy as the dimension 
d of the single system Hilbert space increases. This 
suggests the search for bounds describing the re- 
source limitations in coding more adequately, per- 
haps by taking into account more detailed features 
of the error syndromes than just the maximal num- 
ber of errors. Our construction could be helpful for 
developing and testing such bounds. 

• Non-stabilizer quantum codes can be constructed 
from families of stabilizer codes by taking their 
union p3j ), where one has to require that the pro- 
tected subspaces, corresponding to the codes within 
the family, are mutually orthogonal and that this 
property remains valid after error operations. Ex- 
amples of such non-stabilizer codes are given in 
[ pjj . [25| . In view of our construction scheme, it 
would be desirable to find sufficient conditions for 
a family of graphs such that the union of their cor- 
responding graph codes yields a (possibly more ef- 
ficient) non-stabilizer code. 
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APPENDIX A: 



Proof of Theorem 

tion wp^], defined by 



[II.l 



We first compute the func- 
It is convenient to introduce 



for two subsets K,K' of X U Y the expression 



X r ( 9 \ 9 K ') 



n x( M ) r(fe < fc,) (ai) 

{k,k'}:keK,k'eK' 



where the product is taken over all two elementary sets 
with one element taken from K and the other taken from 
K' . Hence the factor for {k,k'} — {k',k} only occurs 
once within the product. Now we write the integrand in 
(|l2| ) as a product of two terms 



vr\ S x ,g E ,g I }Mh X ,h E ,g I ] 

lX r {h XuE 1 h XuE ) x r (h XUE ,g I ) 



(A2) 



\G\ 



x\- 



X r (9 XuE ,9 XuE ) X V {9 XUE ,9 I ) 



Here only the last factor on the right hand side depends 
on the integration variables g 1 associated with the set 
I = Y \ E. In order to carry out the integral over one 
variable g^, i G /, we select the gi dependent part out of 
X F {h XuE ',/), which is 



{z,i}:z£XUE 



(A3) 



Analogously the gi dependent part of x (g 
the same expression with h replaced by g. Thus the gi 
dependent part of (A2) is 



n x(fc^i)- r(z ' i) x(^,5i) r(z,<) 

{z,i}:zeXUE 



(A4) 



Using the character property of x('i9i) we can simplify 
this to a single factor of the form x(k,gi)- Explicitly, 



X! T v( h o~9j) 
jexuE 



(A5) 



and this sum contains none of the variables associated 
with J, because I (~l (X U E) = 0. The integral over gi 
then gives 5(k), and we find 



|G|i 



XUE uXUE^ 

E UXUE 



w [r ,B] [g"^,h 



x r (g xuE ,g XUE ) 



S(r xuE (h 



XUE 
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XUE 



(A6) 
)) ■ 



Our task is to establish the necessary and sufficient con- 
ditions for this to be of the form 



C(g»,h»)8(h*-g*) 



(A7) 



required by Eq.([13|). 

Now the expression (A6) has the required property of 
vanishing except for g x = h x if and only if this is already 
implied by the vanishing of the (5-function in ( |A6| ), i.e., if 
and only d x ^ implies T XllR {d XuE ) £ . This is the 



first part of the condition in Theorem III. 1 
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From now 



on we assume, 
ix _ 



as we may, that 
implies d x = 0. Then the depen- 
dence of ( |A6| ) on the input variables g x and h x can be 
simplified. The 5-function can be written as 



S(T{ 



xu E (h XuE - 9 XUE )) (A8) 
= 5(T E (h E -g E )) S(h x -g x ), 



x 



because the two expressions are equal for h 
for h x 7^ g x they both vanish by assumption. 
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x 



and 



To simplify the bicharacter quotient in (AC), we use 
Eq. (CO) to write 



X r (h XUE ,h XUE ) (A9) 

= x r (h x ,h x ) x ^h E ,h x ) x ^h E ,h E ). 

With a similar decomposition of x r (<? XuS , g XuE ) we use 



the condition that, wherever the (^-function in (A6) is 
non-zero, we have g x = h x . Hence the X — X factors 
cancel, and we can write the quotient of the X — i?-terms 
as 



X l (h E ,h x ) 

x T {g E M) 



IIx(E r ^ 

jex %eE 



(A10) 



For (|A|) to be of the desired form (|A7|) with C(g E , h E ) 
independent of the X-variables, this expression must 
be independent of all hj,j £ X, whenever S(T E (h E -- 



g E j) 0. But (AlC) is independent of hj if and 



only if '^2 i Tji(hi — gi) = 0. Hence we must have that 
T E (h E - g E ) = implies T% (h E - g E ) = 0. This is 



the second condition from Theorem [II. 1, which we have 



thus shown t o be necessary. Conversely, it is sufficient to 
ensure that ( A10| ) is equal to 1 , and ( A6 ) has the desired 
form (A7) with 



C(g E ,h E ) 



5(T E (h E 



x r (g E ,g E ) 

This concludes the proof. 
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